Traveling waves and geometric scaling at non-zero momentum transfer 
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We extend the search for traveling-wave asymptotic solutions of the non-linear Balitsky-Kovchegov 
(BK) saturation equation to non-forward dipole-target amplitudes. Making use of conformal in- 
variant properties of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) kernel, we exhibit traveling-wave 
solutions in momentum space in the region where the momentum transfer q is smaller than the 
characteristic scale Q of the projectile. We prove geometric scaling in the variable Q/q£l s (Y) where 
Q S {Y) has the same energy dependence as in the forward analysis. 
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(NJ ; I. INTRODUCTION 

Geometric scaling [1| is an interesting phenomcnological feature of high energy deep-inelastic scattering (DIS). It is 
expressed as a scaling property of the virtual photon-proton cross section, namely a 1 (Y, Q) — a 1 (Q/Q S (Y)) . Q is 
the virtuality of the photon, Y the total rapidity and Q S (Y) an increasing function of Y, called the saturation scale 
0- On the theory side of the problem, it is convenient to work within the QCD dipole picture of DIS 3]. . In this 
| framework, the geometric scaling 

> : 
o . 

appears to be a genuine property of the conveniently- normalised dipole-target forward scattering amplitude Af(Y, p) , 
' where p is the size of the dipole and Y its rapidity. In fact, it is known that, if this amplitude verifies the non-linear 
Balitsky-Kovchegov (BK) saturation equation 0, |j| , this property is a mathematical consequence of the high-energy 
behaviour of the solution in terms of traveling waves Q . More precisely, the equation is shown to admit traveling- wave 
solutions which translate directly in terms of the geometric scaling property Q . 

However, the BK equation is written for the dipole-target amplitude in the full 2-dimensional transverse coordinate 
Q_i| plane. Hence, it is supposed to give solutions for the dipole amplitude Af(Y, p, b) where b is the impact-parameter of 
JL . the dipole projectile with respect to the target. Saturation as a function of impact parameter has been studied either 
qj' phenomenological ly B , in the framework of models p|, from a semi-classical approach 0, from numerical studies of 
the BK equation |lOl 1 1 1| or from an analytic point of view [T^. However, referring to solutions of the BK equation, 
the perturbative tail in impact parameter shows up fastly in the BK evolution [T(il | . It appears difficult to be consistent 
in a region where the confinement is expected to dominate (see discussions in Refs.jlj, 

In the present paper we want to tackle the problem of non-forward amplitudes in the saturation regime from a new 
point of view, motivated by (and by extension of) the mathematical properties of non-linear evolution equations, used 
in Refs.jfJ for the forward amplitude. In the present work, we shall stick to the framework of the BK equation but 
our method could be extended to different cases with the same mathematical features. 

In Refs. 6], it was shown that the BK equation for Af(Y, p) can be considered to lie in the same universality class 
than a well-known equation, namely the Fisher or Kolmogorov-Petrovsky-Piscounov (F-KPP) equation Jl5|. It has 
been shown [lq that these equations admit asymptotic traveling-wave solutions whose physical meaning |fj is nothing 
else than the geometric scaling property. 

More generally, if some general features that we shall now point out are realized one can infer the existence 
and the form 6] of traveling wave solutions only from the knowledge of the linear part of the kernel ^^|. Let us 
quote these general features, taking as an example the general structure of the F-KPP equations. It is known that 
they are governed by three types of terms: a "diffusion term" , a "growth term" and a non-linear "damping term" . 
In particular, suppose the "damping term" to be absent, the equation is restricted to its linear part and leads to an 
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FIG. 1: Critical exponent y c for the full BFKL kernel. 



exponential rise of the solution with "time" together with diffusion in "space" . This is indeed characteristic of the 
BFKL kernel which governs the linear part of the BK equation. In more general cases, called "pulled front" cases 
in the literature |17| . the presence of similar ingredients, together with a specific property of initial conditions being 
steep enough to carry along the critical regime of traveling waves, induces the existence and the form of traveling wave 
solutions. The key point is that the main features of these solutions can be determined only from the knowledge of 
the solutions of linear part of the equation, which is an easier task than looking directly for solutions of a non-linear 
equation. 

Our starting point is the knowledge of the exact solutions of the BFKL equation for dipole-dipole scattering, which 
have been obtained using the conformal invariance of the BFKL kernel |2fJ, l2lL |22|. Our aim is the following: 
making use of the powerful mathematical properties of conformal symmetry, allowing to formulate exact solutions 
of the BFKL equation with impact parameter and/or momentum transfer, we can construct the solutions of the 
linear part of the BK equation. We shall then look for the kinematical domain where the ingredients allowing for the 
existence and derivation of traveling waves are present. Namely, schematically (these points shall be discussed more 
completely in the next section): 

• A solution of the linear part of the equation with an exponential growth, 

• A non linear damping effect due to unitarity, 

• A steep enough initial condition. 

By contrast with previous approaches of non-forward amplitudes, we shall not try to get solutions in the full phase 
space (which after all is not necessarily dominated by saturation effects) but shall concentrate on kinematical regions 
where the mathematical requirements for traveling waves are fulfilled. 

The plan of our study is the following. In section [HJ we recall in detail the derivation of the traveling-wave 
speed (i.e. the saturation-scale energy dependence) and of the wave front (i.e. the amplitude) in the case of the 
forward amplitude. In section IIHI we derive the solutions of the non-forward BFKL dipole-dipole amplitudes in 
various representations : full coordinate space, a mixed one where the external particles are coordinate space dipoles 
interacting with a given momentum transfer q, and full momentum space where the external particles have given 
momentum k, Uq. In the next section ITVl we construct the solutions of the linear part of the BK equation and select 
the kinematical domains and representations in which one meets the abovementionned requirements and we derive 
the corresponding properties. Finally, in section we summarise our results and derive their main consequence, 
namely geometric scaling at nonzero transfer. We show how it provides a solution to the puzzling problem related to 
the confinement scale. We point out interesting phenomenological consequences of our solutions. 



II. TRAVELING WAVES IN THE FORWARD CASE 



Let us recall how traveling waves and geometric scaling emerge from the study of the BK equation for forward 
amplitudes. 
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In the large- N c approximation, it is well-known that the high-energy behaviour of the dipole forward scattering 
amplitude off a large (e.g. nuclear) target follows the Balitsky-Kovchegov equation 0,0 where we neglect the impact 
parameter dependence. This equation can be put into the form Q 

d Y Af(Y, k) = a X (~d L )Af(Y, k) - aAf 2 (Y, k), (2) 

where 

X(7) = 2^(1) - V>(7) - - 7) 

is the BFKL kernel and 

Af(Y,k)= ^ ^-J Q [kp)N{Y,p) 
Jo P 

can be interpreted as the density of gluons in momentum space. Here, L — log(fc 2 /fcg) with kg some fixed scale. 
The starting point of Ref.jy] is that, if we expand the BFKL kernel to second order around 7 =§, this equation is 
equivalent, up to a change of variable, to the Fisher-Kolmogorov-Petrovsky-Piscounov (F-KPP) equation [l5| 

dtu(t, x) = d 2 u(t, x) + u(t, x) — u 2 (t, x) (3) 

with t ocY. Looking to the different terms in the right hand side of ©, one can identify a "diffusion term" (d 2 u), an 
"expansion term" (u) and a "damping term" (— u 2 ), as explained in the introduction. 

This equation has been extensively studied and it has been proven that it admits traveling- wave solutions [HI i.e., 
at large time, the solution can be written 

u{t,x) ~ f(x — m(t)) 

t — >-(-oc 

with 

m(t) = vt - tulogf + zt- 1 ' 2 + 0(l/t) , 

where the constants v,w,z can be determined |l7j from the properties of the linear regime. If the initial condition 
behaves like e~@ x with (3 > f3 c — 1, these coefficients acquire critical values, whatever the value of (3 is. In particular, 
the speed v takes the critical value v c = 2. The three coefficients are the only "critical" ones since they do not depend 
on the specific form of the non-linear damping. 

In Ref.@, it has been shown that this kind of feature is expected to be much more general. Let us for instance 
show the properties of the critical regime, when the linear part of the evolution admits a superposition of waves for 
solution. One writes: 

u(t,x)= / £L Uo ( 7)e -T-M7)t = / |X Uo(7)e -7(*w^)+"(7) tj (4) 

J c — zoo Jc-ioo ^27T 

where lo(j) is the Mellin transform of the linear kernel and x w f = x — vt is the position relative to the wavefront. 
Then, the non-linear term drive the solution to the critical behaviour corresponding to traveling waves at large time 

u(t,x) - e -7c*W. (5) 

t — >oo 

The critical exponent 7 C corresponds to a wave having a minimal phase velocity equal the group velocity (see e.g. 

Vc = ^hA = d h){ (6) 

The appearance of traveling waves 10 with critical velocity given by JfjJ only depends on a few general conditions: 

• u = is an unstable fix point of the equation, and u = 1 is a stable fix point. 

• the linearised evolution equation admits solutions of the form It generates a growth of the solution and 
non-linearities damp the solution and saturate it. 

• The initial condition is steep enough. If u$(x) ~ e~ loX , this means that we want 70 > J c - 



FIG. 2: Relevant variables in the collision of two dipoles viewed in the transverse plane. 



In particular, in the case of the 6-independent BK equation, 0^(7) = ax(j) and we find j c ~ 0.6275, as represented 
in Fig^ As noticed in |(| we expect 70 = 1 due to QCD colour transparency and, therefore, the QCD evolution in 
rapidity will asymptotically reach the traveling wave regime (jSJ) which in turn corresponds to geometric scaling 

where the saturation scale Q 2 S {Y) grows as k^e VaY and the critical speed v c is 4.8834a. 

In this paper we investigate in which kinematical regions one has similar properties from the BK equation in full 
phase-space. We shall first discuss under which circumstances the BFKL dynamics leads to solution of the form 
(0J. We shall then show how we can obtain solutions for the linear BK equation and how non linearities lead to the 
formation of traveling waves and to geometric scaling. 



III. LINEAR BFKL DYNAMICS AT NONZERO MOMENTUM TRANSFER 

If we now take into account the 6-dependence in the BK equation, we have to study the asymptotic solutions of 

cVAA(x, y) = [ d 2 z (x ~ y)2 [AT(x, z) + Af(z, y) - A^(x, y) - Af( X , z)Af(z, y)] , (8) 

where jV(x, y) is the conveniently-normalised dipole-target scattering amplitude, x and y are then the transverse 
space coordinates of the quark and antiquark constituting the dipole. This equation does not depend explicitly of the 
target whose centre of mass defines the origin. This non-linear equation corresponds to resumming QCD fan diagrams 
in the leading- logarithmic approximation pg. 

Our main tool is to use the knowledge of the exact BFKL solutions 0, 120I l2lj for dipole-dipole scattering. We 
shall study the amplitude f(pi, p 2 , Py , P2') ( see Fig|2J) where we identify p x and p 2 with x and y . We impose 
that one dipole is bigger than the other (p = \p 1 — p 2 \ <C po = \py — Pa' I ^ Indeed, if we want to 

use the dipole-dipole amplitude in order to study the dipole-target amplitude within the BK framework, we need to 
extract solutions of the linear part of the BK equation |JS} from the solutions of the BFKL equation. This is done 
by integrating out the target impact factor, and gives relevant results provided the convolution with the large dipole 
factorises from the smaller one. This separation criterium closely depends on the conformal invariance properties of 
the BFKL equation which is one of the basic properties of the non-forward linear BFKL dynamics |l9j| . 

Since our discussion may lead to different theoretical and phenomenological features, we shall perform our analysis 
in different spaces: the coordinate space, the mixed representation where we keep the dipole sizes in coordinate space 
but use the momentum transfer q instead of the impact parameter, and the full momentum space. 

A. Coordinate space 

Let us start from the well-known BFKL amplitude 0] 
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where s = e Y is the centre-of-mass energy squared and q is the momentum transfer. One writes |19| 

z/ 2 ^ i 

2tt(^ 2 + 1/4) 2 u - ax{v) 



f u (Px,P2,Pv,Pv) = I ^,,1 7Tuv> - — \tt^ I <PrE v {p x -r, P2 -r)E v (p v -r,p 2 ,-r) CM 



where 



^(pi,P 2 ) 



, l+2iy 

iPl - Pl\ 



\p\\m\ 



are the conformal eigenfunctions of the BFKL kernel with vanishing conformal spin quantum number (n — in the 
usual terminology 19]). It corresponds to the dominant contribution at high energy. $p (resp. $t) are the impact 
factors describing the coupling of the BFKL pomeron to the projectile (resp. target), depending on the momentum 
transfer q and on the dipole sizes p = p x — p 2 and p$ = p v — p 2 ,. 

The Dirac distribution which appears in the definition of / w ((j 2 ) comes from the fact that, the global system being 
invariant under translation, the expression in the right-hand-side does not depend on Px+Pz+Pv +P 2 '- To be more 
precise, the relevant variables of this problem are the sizes of the two dipoles, p and p , and the impact parameter 
b = ^(pi + p 2 ~ Pv ~ Pt\ as shown on Fig|21 Combining these expressions and using 1 7 = i + iv instead of v, we 
obtain the following expression for the amplitude 

A ^1 2 ) = TTvJ ^Wpoe 4q b $p(P,q)<i>T(Po,q) / 2+IO ° ^e a *(^r(p, Po ,b), (10) 
where we have introduced the function / 7 as the result of the integration over r in 

P(p, Po ,b)= J A\E\p 1 -v,p 2 -x)E\ Pv -v,p 2 ,-v). 

Using the complex representation for the two-dimensional vectors 2 (in this section, we keep bold characters for vectors 
and ordinary ones for complex numbers), it has been shown pol l22j ] that the result of this integration is 

rip, Po, b) = 2n ° 7 2 . |z| 27 2 F 1 (7, 7; 27; *)aFi(7, 75 2 7 ; z) + (7 <-» 1 - 7), 

where 2F1 is the Gauss hypergeometric function [25|. the pre- factor 

r(i-7)r (7- ±) 

and 

P12PV2' _ 4pp 



Pll'P-22' 46 2 - (p - po) 2 



(12) 



is the anharmonic ratio. This is a remarkable property due to conformal invariance that the BFKL pomeron exchange 
in coordinate space depends on p, po an d b only through the anharmonic ratio z. 

Since we are interested in the situation where one small dipole of size p scatters on a larger one of size po, i.e. 
|p| <C |po|, we can expand / 7 (p, Po,b) in series of p/po- In that limit, the hypergeometric function goes to 1 and we 
obtain 



/ 7 (p,p ,6) 



7 2 (1"7) 5 



4ppo 



27 



+ (7^1-7). (13) 



1 In the following, we shall use indifferently 7 or v as superscripts, e.g. E v and E 1 denote the same function. 

2 The complex representation of the vector x = (xi,X2) is 

x = xi + iX2, x = x\ — iX2. 
In this section, we shall explicitly use \x\ when the modulus of the vector has to be considered. 
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Before going any further, one has to point out that the corrections to this expression are of order (p/po) 2 but reduce 
to (p/po) 4 if we integrate over the phase (f> of p (see FigEJfor the kinematics). 

The expression still depends on the angle 4>q between po and b. Since this is not relevant for phenomenological 
studies, it is interesting to integrate this result over the phases <f> and <f>o of p and po (one may assume that b is real). 
One obtains 



(/ 7 (p,Po,&))«^ 



4|p||p | 



2", 



7 2(1_ 7 )2 



|po| 2 -4|6| 2 



7 -l 



I P0| 



i6i&r 



I Pol 



16|6| z 



+ (7^1-7), 



where P^-i(x) is the Legendre function of the first kind [25| . 

Finally, if we additionally take b going to zero in this result, we obtain 



1 



7 2(1_ 7 )2 

which exhibits the same power behaviour as in the forward case. 





Ap_ 


27 


ip 


2-2 7 \ 


( C 7 




+ ci_ 7 






Po 


Po 





(14) 



(15) 



B. Mixed space 

It is useful to introduce a mixed space representation |19| in terms of p, po and q. One obtains from i|l(J|) 



A(s,q 2 ) 



(2k) 



d 2 pd 2 p Q ^ P {p,q)<S? T ( Po ,q) j —f»(p,p ) 



dv 



where 



f:(p,p ) = 



d 2 b 



Jq.h pv 



f v (p,P ,b) 



IPI iPol 



(2tt) 2 j ; 16(^ 2 + 1/4) 2 '' 

acquires a factorised form with 
Using again the complex representation for the two-dimensional vectors, one finds (see |20|') 



(16) 



(17) 



E»{p) = |<?r 2^'T 2 (1 - iv) [j- iv (f ) J- iv 



— J, 



m (T) Jiv (T) 



and 



E"(p ) = E- v (p ), 



where J a (z) is the Bessel function of the first kind. Using 7 
into a (7 <-» 1 — 7) symmetry and we finally get 



iv instead of v, the (y «-> — v) symmetry turns 



fq(p,Po) 



lPPo| 



16 7 2 (1 - 7 ) 2 



7 r 2 



(f)^(f) 



(?) 



(18) 



"3 V 4 / 



•2 V 4 / 



V 4 / 



V 4 / 



Let us notice that this expression has the crucial feature that it is factorised in p and po by contrast with the b 
representation. This property, inherent to the fact that we use the momentum transfer q instead of the impact 
parameter 6, will be of prime importance for obtaining the solutions of the BK equation. 

As explained before, we are interested in the situation where one dipole is much larger than the other one (po 3> p). 
In order to expand the Bessel functions in series of their respective arguments, we shall consider three situations: 
> p > 1/<Z, Po S> 1/q 3> p and 1/q po p. The results then follow from the asymptotic expansion of the 
Bessel functions 



0) J„ (z) - J_ M (z) J_p (z) 



1 



r2(i + M ) 



2// 



1 



r 2 (i-/i) 



-2/i 



7T Z 



sin(^7r) cos [29?e(z)] . 



(19) 
(20) 



Using these expressions we obtain the following results for the three relevant cases. 



7 



Case 1 : l/q > p > p. 
Using (P)l and flg) , we find 



fq(p,Po) 



16 7 2(l- 7 )2 r2 (2 7 ) r "(2 +7 



1 



r 2 (l 



7) 



1-27 



1 



r 2 (^ + 7) 



P5 



27-1 



r 2 (H7) 



27-1 



r»(f 



7) 



1-27 



In particular, if we take q — > in this expression, we recover 

q^o \pp \ 



fq{p, Pa) 



16 7 2 (l- 7 ) 2 





1-27 

+ 




2 7 -l\ 


p_ 


p_ 


Po 




Po 





(21) 



• Case 2 : p > l/q > p. 
Due to the asymptotic expansion l|2(J|) . the po part will not depend on 7 anymore and we obtain 



fa(P>Po) 



— 1 C0s( 7 7r) 
"2V 7 2(1_ 7 )2 



"Ke(p g)" 


1 


KH 


pq 




P<7 


2-27" 


2 




8 


8 





(22) 



Note that this amplitude still depends on the angle i[iq between q and po. If we integrate out this angle together 
with the dependence on the angle ip between q and p, we have 



</ 7 (P,Po)) 



- C0S( 7 7I") 



^,V>o 7 2 (1 — 7 ) 2 



\P00\ _ 7T 

2 4 



1 



7T 3 ipo? 3 ! 



r-'l|-7 



2-27' 



(23) 



• Case 3 : p a ~> p~> l/q- 
In this case, both the p- and po-dependent parts of (|18|l involve expansion (|20|l which gives 



— 4 COS 2 ( 7 7r) 



1 



7r 2 7 2( 1 _ 7 )2^ \ 2 

and, averaging over the angles ip and ^>o> 

1 



r* U+7 ^ --7 t^cos 



1 



5Re(w) 



3?e(/9 g) 



</ 7 (p,Po)) 



-16 cos 2 ( 7 tt) 2 
7T 4 7 2 (1 - 7 ) 2 \ 2 



7 P 



cos 



IP9I 
2 



cos 



(24) 



(25) 



C. Momentum space 

The situation in momentum space is very similar to the one in mixed space. If we still neglect the contribution of 
higher conformal spins, we find 



A( S ,q 2 ) = is J d 2 fcd 2 fc $ P (k,q)$ T (k ,q) J £ e & ^ Y f(k, k , q), 



where <J>p and <&t are the impact factors in momentum space 

$ P (k,q)= J ^ e -(-c/2).P $p(p)q 

and / w (k, ko, q) is the Fourier transform of f" (p, p ): 

r(k,k ,q) 



(2tt) 4 (^ 2 + l/4) 2 



(26) 



(27) 
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The detailed calculation of this integral is presented in Appendix ^] and the relation with conformal invariance is 
developed in Appendix^] The final expression, using again the complex representation (see footnotc|2l and replacing 
v by 7 = h + iv, is 



P{k,k ,q) 



-1 



sin 2 (771") 



(2tt) 4 16 7 2 (l-7) 2 

r 2 (i + 7 ) a 2-/- 1 



7 r 1 



\kkr, 



r 2 (^ + 7) 

r 2 (i + 7 ) 
r 2 (^ + 7) 



JL 

4k 



2 Fi 7, 1 + 7; 2 7 ; f 2 F X 7, 1 + 7; 2 7 ; i - (7 1 - 7) 



4fc n 



27-1 



2F1 (7, 1 + 7; 2 7 ; |j 2F1 f 7, f + 7; 2 7 ; ^ 



(7^1-7) 



(28) 



As in the mixed space representation, this expression has the remarkable property of being factorised in k and ko. 
This appears from the fact that we use the momentum transfer q and is lost if we go back to impact-parameter b. 

We can now consider these expressions in the case where one dipole is much harder than the other, corresponding 
to k 3> ko- Thanks to the fact that p is factorised in q/k and q/ko, we again need to consider three situations: 
k 3> ko >• q, k ^> q ^> ko and g > <; > fco- When q is small compared to k, it is sufficient to consider 



2F1 (7, 1+7; 2 7 ; |) 



q\ |fc|>kl 



1 



The behaviour at large q is more complicated. We show in Appendix [U] that, averaging over the angle 9 between k 
and q, we have 



/r 2 (i- 


-7) 


q 


27 


V 2 {\- 


f 7) 


4k 





2F1 (7, 1 + 7; 27; f ) 2F1 (7, 1 + 7; 27; I) - (7 <-> 1 - 7) 



l fc l«H d 2(1 7)2 cos(7^) 
sin(77r) 



log 



(29) 



Before considering the three cases in more details, one may note that the q~ 3 behaviour could have been anticipated 
from the fact that, in that limit, the Fourier transform of E^(p) does not depend on k. Note however the logarithmic 
term arising from the hypergeometric function. 

Let us now consider the results for the three different cases. 

• Case 1 : k 3> k » q. 
In this case, both parts of the expression show a power behaviour, 

-1 sin 2 (77r) 



P(k,k ,q) 



7T 4 4 7 2 (1 - 7 ) 2 



7 I" 



1 



+ 7 



r 2 (i + 7 ) q_ 2 7-i T 2 (2- 7 ) 
4k 



r 2 {\ + 7) 



(27T) 2 |fcfc |^ 



r 2 (f-7 

1 27 sin 2 (777) 



Q_ 

4k 



1 

\kk 
1-27 



r 2 (i + 7) 
r 2 (i + 7 ) j/m, 



r 4 (i + 7 )r 2 (|- 7 ) 



27-1 



r 2 (2-7) 
r 2 (I - 7) 



(30) 

1-27 



16tt 2 7 2 (1 - 7 ) 5 



r 2 G + 7) 



l6kko 



27-1 



(7 



1 



4fc 

-7) 



In the distributed product, we clearly see that, taking the limit q — > 0, we recover an expression depending only 
on the ratio between the momenta of the two dipoles 



/ 7 (Mo,<z) 



1 



1 



(2tt) 2 |fcfc r 



k 


1-27 


k 


2 7 -l\ 




+ 






k 


k 





(31) 



• Case 2 : k ^> q 3> k . 

At intermediate values of q, the dependence on ko goes away and we recover an analytic power behaviour in 7 
which, as we shall see in the next section, will lead to traveling waves and geometric scaling 



P(k,k ,q) = ^sincarr)!*!---. )r-'(- 



log 



\kq\ C 



T 2 (l+7) 

r 2 (l + 7) 



4k 



27-1 r 2 (2- 7 ) 



r 2 (l"7) 



q_ 

4k 



1-27 



(32) 
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Case 3 : q >• k » k . 
This limit gives 



P(k, ko, q) = -£ cos 2 ( 7 ^)r 2 (~ - 7) T 2 (\ + 7 ) 7 2 (1 - 7 ) 2 ^ log 



log 



(33) 



IV. TRAVELING WAVES AT NONZERO TRANSFER 



In this section, we investigate the consequences of our formulae for linear BFKL dynamics on asymptotic solutions 
when the non-linear effects of the BK equation are switched on. The key point is that linear BFKL dynamics can 
provide solutions for the linear part of the BK equation. Indeed, when factorisation between the target and the 
projectile is fulfilled, the impact factor of the target can be integrated out keeping a linear evolution governed by the 
same kernel. It is thus a solution of the linear part of the BK equation which should depend only on the kinematic 
variables of the projectile. 

In the next step, the asymptotic solutions of the BK equation can be deduced from the knowledge of its linear part, 
using quite general arguments which were recalled in Sections [I] and |nj For each of the kinematical configurations 
introduced in the previous section, we check whether the conditions for the existence of traveling waves are fulfilled 
and, in these cases, derive their expression. 

Let us first concentrate on the full momentum representation J2HJ)- We define the function 



/(k,q) = J d 2 fc $ T (k ,q) J |Le a ^r(k,k ,q), 



(34) 



depending only on k and q. It is remarkable that in this representation, the exact factorisation property of / 7 (k, ko, q) 
(see (|23) gives rise to the simple expression 



where 



/(k,q) = J |le a ^ y ^(q)r(k,q), 
q)= / ^fco^ e l(k °" q/2) - P «^^(Po)'i>T(ko,q) 



(35) 



factorises out the target dependence, defining the initial condition, and 



2 4-6 7 r 2(3 



V ■ / s\Q 
sin( 7 7r) — 



27-1 



r 2 (i + 7 ) 
r^i + 7) 



4k 



27-1 / q 
2FW7 + 1, 7 ; 2 7 " 



! } 2 Fi(7 + l,7;27;| 



(7 



(36) 
7) 



is calculated in Appendix iBl and enters in formula 128(1 . It is obvious from the structure of l|34|) . a linear superposition 
of eigenfunctions of the BFKL kernel, that it provides a natural basis for the solutions of the linear part of the BK 
equation in momentum space. 

Therefore, we look for the kinematical domains where / 7 (k, q) takes the appropriate exponential behaviour, cor- 
responding to the wave superposition property discussed in section [n] It will lead to traveling- wave solutions for 
the full BK equation. 

We observe that, both in the cases 1 and 2 (formulae (|31|l and 132f) '). i.e. in the limit k^q, the solution of the 
linear part of the BK equation can be recast, using the symmetry 7 «-> 1 — 7, under the form 



where 



L = log 



expresses the leading power-like behaviour at large k/q and all remaining factors in 7 and q has been reabsorbed in 
4>"'(q). This clearly emphasises the fact that the region of interest is k^>q, independently of kg. By contrast, in the 
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case 3 treating the large-g expression, we see that equation (|33ll . though being factorised, does not fulfil the same 
property. We thus do not expect traveling- wave solutions in this limit. 

Let us now include the effect of the nonlinear term in the BK equation. Due to colour transparency (/(k, q) ~ k~ 2 
for k ^S> q), one can apply the general arguments exposed in Section [H] Therefore, the solution of the BK equation 
reaches asymptotically the traveling-wave structure and exhibits geometric scaling. The saturation scale 



Qi(Y) = q ^ s (Y) 



(37) 



q 2 exp 



^r-A log(F) 

7c 27 c 



has the same rapidity dependence as in the forward case but is proportional to the momentum transfer. The critical 
exponent j c is still given by equation ©. The asymptotic form of the amplitude, i.e. the wavefront, can be written 
at small q/k 



/(k,q) 



Ta ^ fa) lo § 



q 2 n 2 {Y) 



k 2 



q 2 n 2 (Y) 



exp 



1 



2a X "{lc)Y 



log 2 



fc 2 



q 2 n 2 (Y) 



(38) 



It is interesting to see how we recover the forward limit q — > recalled in Section^] This can be done remarking that 
formula H30|) . obtained for BFKL dynamics in the case 1 (k ^> ko ^> q), introduces and additional factor (q/ko)~ 21 . 
By recombination of the q/k and q/kg factors, a different factorisation appear where fco substitute to q as the reference 
scale, as clearly seen in (|31ll . Inserting (|31|) in (|34p. by straightforward algebra, it is easy to see that the result can 
be cast under the form 



/(k,q^0) 



1 

k 2 



dry 
2i^ 



f (Q T )e 



where 



u = log (m 

and Qt is a scale typical for the target, defining the initial condition for the traveling- wave solution. This corresponds 
to the solutions of equation J2J. 

This discussion can be translated to the amplitude in the mixed space representation of section IIII Al defining 
similarly 



/(P,q) 



'Po 



{2-Kf 



*T(p ,q) 



d'y 



fq(P,P ), 



where f^(p,p ), given by l|17|l . is the solution of the BFKL equation in the mixed representation. As for the case 
of the momentum space, the factorisation of the target and projectile dependences in the BFKL amplitude is the 
important property, leading to traveling waves when we include the nonlinear effects in the limit p <C 1/q. 

The situation is not as simple in coordinate space. Indeed, it is hard to introduce an amplitude depending only on 
p and b from equation l(TT)|) . This is closely related to the fact that / 7 (p, p ,b), depending only on the anharmonic 
ratio, cannot be factorised. In the limit where the BK equation is valid (po ^§>b, p), its linear part gives solutions of 
the form (see (EEl) : 



/(P,b) = J 



2a, 



2m 7 2 (l- 7 ) 2 



, x i,n ( 

Po 



27 



where po is the typical size of the target. As explained in Section [HJ this leads to traveling- wave solutions for the 
BK equation J\f(p,Y) —Af(p Q S (Y)), where all trace of the impact parameter has disappeared. Hence, it gives no 
information upon the b dependence. We could instead consider equation (|14fl when the scaling variable ppo/\Po ~ 46 2 | 
is large, i.e. p ~ 2b. However, it is not clear whether the BK equation makes sense physically in that limit where the 
dipole hits the target in its peripheral region (po ~ 2b). As explained above, these difficulties arise from the fact that 
p and po are mixed non-trivially in the anharmonic ratio (|12f) . Hence, an interpretation in term of the BK equation 
is problematic. By contrast, the situation in mixed and momentum spaces does not suffer this inconvenient, due to 
the nice factorisation property in p and po, or in k and fco- 
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V. CONCLUSIONS AND PERSPECTIVES 



Let us first emphasise the main results of this paper. We show that traveling- wave solutions of the forward BK 
equation can be extended to the full equation including nonzero transfer q, provided that k 3> q, where k represents 
the scale of the projectile. The saturation scale lj3T|l has the same energy dependence than in the forward case but 
is now proportional to q. When this scale is large w.r.t. the scale, e.g. Qt, characterising the target, the saturation 
scale thus becomes independent of the details of the target. But, when the momentum transfer q is small w.r.t. the 
target scale, Qt substitutes to q in the saturation scale as expected from the forward case analysis. Our results 
show that non- forward BKFL evolution with momemtum transfer q in a presence of saturation at the scale Q S (Y) is 
equivalent to forward BKFL evolution in the presence of saturation at the scale qil s (Y)/QT- Note that the similarity 
of a momentum transfer with an absorptive boundary in BFKL evolution has been noticed in Ref.|2J|. 

As a consequence of the existence of traveling- wave solutions, we derive the asymptotic form of the solution, i.e. 
the wave front (|38f) . at small q/k. It appears as an expression of the same form as in the forward case with two 
modifications: the saturation scale is now the proportional to q and the pre-factor, related to the initial conditions, 
is now q-dependent. 

From these considerations, we can predict that the BK equation implies the extension of geometric scaling at nonzero 
momentum transfer. Indeed, using formula l|26|) and noting that, at large k, the impact factor of the projectile is 
expected to scale with the ratio k/Q where Q is a typical hard scale of the projectile. We see that the amplitude 
satisfies the geometric scaling under the form 



A(s,q 2 ,Q 2 ) ex is^(q)f I ^ \ . (39) 

All these results apply also in the mixed-space representation depending on the projectile size p and the momentum 
transfer q. However, we find that these properties seem not to be easily expressed in terms of the impact parameter. 

Technically, the method used consists in three main steps. Firstly, we analyse the form of the BFKL solutions in full 
phase space and in different pertinent limits. Secondly, we use conformal-invariant properties of the BFKL dynamics 
to look for factorised solutions which then lead to solutions of the linear part of the BK equation. Thirdly, we infer 
the traveling-wave solutions, induced by nonlinear effects, in the kinematical regions and for initial conditions where 
the universality properties of the BK equation apply. 

Some comments are in order. Our formula i|38[l s ugg ests a solution to the puzzling problem of the compatibility 
of the BK equation with the confinement scale 0, llCt ITU IlsI IT^ | . The key point is to address the problem 
in momentum space where the non-perturbative dependence is factorised as clearly seen from the factor (jp" (q) in 
equation H38|) . which may be characterised by a confinement scale. Fourier-transforming this result back to impact 
parameter space will break this factorisation property. 

The properties of the BK solutions at nonzero momentum transfer suggest to analyse the high-energy behaviour of 
suitable experimental processes. The electroproduction of p-mesons at nonzero transfer has been already analysed in 
impact parameter space 7j and our solutions of the BK equation incite to perform the analysis in momentum space. 
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APPENDIX A: CALCULATION OF /" (k,k ,q) 
In this appendix, we compute the Fourier transform of 

W = ^ M" 2 " r 2 d + „) [j, (f ) Jn (f ) - J- M (f ) J_, (f )] , (Ai) 

with fi — —iv. It is of course sufficient to compute the following integral 



I 
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where v = \k — q/2\ and <j> (resp. if) is the angle between k — q/2 and p (resp. q). This integral has been regularised 
at infinity by introducing a factor p a with —2 < a < —3/2. If we expand the Bessel functions in series and use 

d e i[p«cos(«-m0] = 2 7 re" rm/2 J m (H, 



we find 



j!fc!r(i + M + j)r(i + M + k) J Q 



j,k=0 



dpp 2 +a+ 2» + 2 j+2k J 2{ ,_ k)ipv) . 



The integration over p can be performed using 

poo i /2\ 13 



/3 + m\ f ^ — m 



j3 — m 



This leads to a factorisation of the integral in conformal blocs: 

" 'i + a + 2p^ 



I = sin 



f(k,q)f(k,q), 



where, using the doubling formula, 



3=0 



j\ r(i + M + i) 



i4, \ ^ +2j 



f + f +M+2J" 



4k-2qJ \\k-q/2\ 



r (l + f + M) 

r(i + /x) ■ 



:Fi 



3 a p 5 a p 

4 + 4 + 2'4 + 4 + 2' 



1+W 



2fc - g, 



The hypergeometric function can be simplified (see e.g. Ref.[25|]), and, taking a to zero, we finally obtain 



I = — cos(/x7r) 



r 2 (l + ^) 
r 2 (i + /i) \\k\ 



Ak 



2/i 



2-T 1 



3i „ g 

- + p, - + p: 1 + 2p; — 



2 Fi [ l+p ,\+p;l + 2p-l 



The final result follows from adding to / the corresponding expression with p — > —p and using i|Alfl . 



(A2) 



APPENDIX B: FOURIER TRANSFORM OF E v q {p) 

In this appendix, we show how the Fourier transform of E^(p) can be calculated using conformal invariance. That 
will exhibit the link between /^(k, kp,q) and matrix elements of SX(2,C) j2||. One has to deal with 



J Spe^^E^p) = ^J d 2 Pl d 2 p 2 e ^V k -^V - p 2 r\ Pl \-^\p 2 \ 



-27 



where c 7 is given by l|ll|) . 

To compute this integral, one can switch to complex coordinates and perform the following changes of variables: 
first P2^u = pij pi, then p\ — > w = pi(k+u(q— k)), and finally u— *v — u/(u~ 1). One obtains 



Then, using 



/ d 2 w e me ^\w\ 2 - 2 ^ = ^2 4 - 2 ^{ 2 7 j 

7 r (7-i) 
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and the following result [27| 

' Sv M 2Q1 - 2 |1 - „|2&i-2«i-2|i _ vx \-*>° = n 



r(oi)r(6i-oi)r(i-bi) 
r(i-oi)r(i-6i+ai)r(6i) 



2 Fi(a , a x ; 61; a;) 2 Fi(a , 01; &i; x)+ 



■ ^o^rll^ 6l ^ a0) N 2 ' 2bl 2J 1 iK-bi + l,«i-bi + l;2-6 1 ;x) 2 F 1 (a -& 1 + l,a 1 -& 1 + l;2-b 1 ;x) 
r(2-6i)r(a )r(6i-a ) 



one easily recovers 
d 2 p e 

"r 2 f# 



'»^^EJ(p) = I M 2 - 2 3 - 6 -r 2 (l - ft,) ooB^ji 



r 2 (i + w) 



4fc 



,3 1 q\ (3 1 g, 

2F1 - + iv, -+iv;l + 2iv\ - 2 Fi - + iv, - + iv; 1 + 2ij/; •=• ] — (1/ — » -1/) 



The same integral appears in the calculation of matrix elements of 5X(2, C) |2fi| and more generally in conformal field 
theory [27|- 



APPENDIX C: ASYMPTOTIC EXPANSION AT LARGE q 



We want to emphasise the main steps yielding the asymptotic behaviour (|29|l . We shall start by the asymptotic 
expansion of the hypergeometric function at infinity {z = q/k ^> 1) 

2 2 7-i r (- + y) ( \os(-z) ( 1 

2 F 1 (7 > 7+l;2 7 ;z) = - } 2 , 1 {-z)^ 1 - 7(1 - 7 )-iL^ 2Fl 7 + 1, 2 - 7; 2; - 



A r(i + 7 ) 



fe=0 



k\(k + l)\ 



When combining the z and z parts, most of the pre-factors cancel: 



r 2 (i + 7) 
r 2 (l + 7) 



27-1 



2 F 1 ( 7 ,7 + 1; 2 7 ; z) 2 F 1 ( 7 , 7 + 1; 2 7 ; z) = —- [S 7 (z) + A^z)} [S 7 (f) + A 7 (z)] 

n\z\ 



where 



z 

A y (z) = \ | 7 (1 - 7) [^(7 + 1) + ^(7 - 1)] + (7 + 1)7(7 - 1)(7 ~ 2) Hl + 2)+ z Hl 2) + 0(1/ z 2 ] 

In this splitting we have put in S-y the terms which are invariant under the replacement 7 — ► 1 — 7 and the remaining 
ones in A™. Therefore, if we consider the full ^-dependent term, we easily get 



r 2 (i + 7) 
r 2 (^ + 7) 



27-1 



2Fi(7,7 + 1; 27;z) 2 Fi(7,7 + 1;27; z) - (7 -> 1 - 7) 



= 4n fa(z) [A^z) - A 1 _ 1 {z)] + S^z) [A,(z) - Ar_ 7 (z)] + \A,(z)\ 2 - |^- 7 (z)| 2 } . 

Obviously, the leading term in that development is proportional to z^^z -1 , i.e. to ~ + |, which vanishes if we average 
over the angle between the two vectors. With less restrictions, one can also remark that this term is antisymmetric 
under the replacement k — > —k. 

We thus go to the next order which finally gives a leading contribution of the form 



1 



7 2 (1 - 7 ) 2 [V>(2 - 7) + ^(-7) " <K7 + 1) " ^(7 " 1) 



log Hz 



2 7 2 (l- 7 ) ; 



cos(77r) 1q S [\ z 
sin(77r) \ z \ 3 
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